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We introduce time-dependent, generalized factorial cumulants C™(t) of the full counting statistics 
of electron transfer as a tool to detect interactions in nanostructures. The violation of the sign 
criterion (—1 ) rn ~ 1 C™(t) > 0 for any time t, order ra, and parameter s proves the presence of 
interactions. For given system parameters, there is a minimal time span t m i n and a minimal order 
m to observe the violation of the sign criterion. We demonstrate that generalized factorial cumulants 
are more sensitive to interactions than ordinary ones and can detect interactions even in regimes 
where ordinary factorial cumulants fail. We illustrate our findings with the example of a quantum 
dot tunnel coupled to electronic reservoirs either in or out of equilibrium. 
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I. INTRODUCTION 

The stochastic nature of electron transfer in meso¬ 
scopic conductors gives rise to both thermal and shot 
noise. 1 A deviation of the shot-noise power from the 
value expected for Poissonian processes of uncorrelated 
charge transfer indicates correlations. In particular, an 
enhanced, super-Poissonian shot-noise power is a clear 
signature of the presence of interactions. 2,3 Various sce¬ 
narios for super-Poissonian shot noise of currents sus¬ 
tained by electron tunneling have been studied both the¬ 
oretically 4-10 and experimentally. 11 17 

Full information about the counting statistics of charge 
transfer 18,19 is contained in the probability distribution 
P/v(t) that N charges have been transferred through 
the system in time t. From its Laplace transform, the 
moment-generating function A4(z,t) := e Nz 
one can derive moments M ni {t) := t)\ z= o and 

cumulants C m (t) := <9™ln M(z, t)\ z= o as the rath deriva¬ 
tive (ra > 1) with respect to z taken at z *= 0. The 
rath moment is the expectation value of the rath power, 
M m (t) = ( N m ) (t), with (• • •) := Eat ''' p N(t). For at 
least two reasons, however, it is advantageous to study 
cumulants C m (t ) = ((7V m ))(£) instead of moments. 20 
First, if the charge transport can be separated into sta¬ 
tistically independent subprocesses, the cumulants of the 
total transferred charge are simply given by the sum 
of the cumulants of all channels. Second, in the long¬ 
time limit, all the cumulants C m (t ) oc t grow linearly 
in t while the moments M m (t) oc t rn grow with differ¬ 
ent powers. 19 Higher-order cumulants have been calcu¬ 
lated for electron transport through various interacting 
systems. 21 38 Experimentally, the cumulants of electron 
tunneling through quantum dots have been measured up 
to 20th order by monitoring the charge occupancy of the 
quantum dot via a capacitively coupled quantum point 
contact. 39-47 

An enhancement of the shot-noise power due to in¬ 
teractions can be quantified by the Fano factor, i.e., 


the ratio of the second to the first cumulant. This 
raises the question whether there is additional infor¬ 
mation about interactions contained in the higher-order 
cumulants ra > 3. Higher-order cumulants have been 
shown 43 to oscillate, as a universal feature, as a func¬ 
tion of any system parameter, measurement time £, or 
order ra, independent of whether or not interactions are 
present. As an alternative probe of interactions, factorial 
cumulants have been suggested instead. 48,49 The latter 
are defined as cumulants (( N ( m ))) of the factorial power 
N := N(N - 1) • • • {N — m + 1). For noninteract¬ 
ing systems, factorial cumulants do not change sign as 
a function of any system parameter or time t and dis¬ 
play only trivial sign changes (—l) m_1 as a function of 
ra. Any deviation from this behavior proves the presence 
of interactions. Interactions are, however, only a nec¬ 
essary but not a sufficient criterion to obtain nontrivial 
sign changes of factorial cumulants. In fact, the results 
of a recent measurement of hole transfer through an in¬ 
teracting quantum dot 47 could be explained within an 
effectively noninteracting model. 

II. GENERALIZED FACTORIAL CUMULANTS 

In this paper, we present a more sensitive and versa¬ 
tile indicator of interactions based on generalized facto¬ 
rial cumulants. To this end, we define the generalized 
factorial-moment generating function 

oo 

E {z + s) N P N (t) 

Ms (z,t):=™- -, (1) 

E s"P N (t) 

N =0 

where N is the number of transferred charges, count¬ 
ing the charges that, say, leave the central part of 
the system (e.g., quantum dot) into some leads while 
not counting those entering; therefore, by definition, 
N > 0. From the generating function, we obtain gen- 
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eralized factorial moments M™(t) := d™M s (z,t) \ z= q 
and cumulants C™{t) := <9IJTn M s (z, t)\ z= Q. Defin¬ 
ing an 5-dependent expectation value (• • •) s (t) := 

• • • s N P N (t)/ s N P n (£), one can show that 

M™(t) = (s- m N^) a (t) and C s m (t) = {{s- m N^)) s (t), 
where (( N rn )) s is related to ( N rn ) s via the same recur¬ 
sive relation that relates ((TV m )) to (N m ). 20,50 The fac¬ 
torial cumulants are recovered by setting 5 = 1, i.e., 
Cf (t) = ((A^M))(t). 

The generalized factorial cumulants with real s can be 
expressed in terms of the zeros Zj(t ) of A4 s (z,t) with 
degeneracies ay via the formula 51 




otj cos [m arg(— Zj)\ 


( 2 ) 


The zeros are either real or appear in complex conjugated 
pairs. By varying 5, their positions are simply shifted in 
the complex plane, i.e., Zj(s) + 5 is independent of 5 . 

For a noninteracting fermionic system, the generating 
function has the form M s (z,t) = (][[[■ [1 ~ Pj PPj( z + 
s)])/E N s N P N (t), where pj with 0 < pj < 1 is the prob¬ 
ability for a single-particle transfer to occur. 52,53 Hence, 
all zeros Zj = 1 — 1 /pj —s lie on the real axis with — Zj > s 
and positive ay. If 5 > 0, then arg(— Zj) = 0, which fixes 
the sign of all generalized factorial cumulants, 

(_l)m-! C m( t ) >q. ( 3 ) 


If 5 < 0, this criterion still holds for all even orders m. 
Any violation of this behavior is a clear indication of 
the presence of interactions. Equation (3), which holds 
for any noninteracting fermionic system constitutes the 
main result of our paper. For 5 = 1 and m = 2, it reduces 
to the well-known criterion that a super-Poissonian Fano 
factor, C 2 (t)/C 1 (t) > 1, indicates interactions. 

What are the conditions on £, m, and s to observe a 
violation of Eq. (3) and thus prove the presence of in¬ 
teractions for a given system? First, at very short times 
£, all charge-transfer processes become independent from 
each other. Thus, all zeros of M s (z,t) are on the real 
axis and Eq. (3) holds. As a consequence, t must be 
larger than a minimal time span £ m i n . Second, higher- 
order generalized factorial cumulants are favorable since 
in Eq. (2) deviations of arg(— zj) from 0 are amplified 
by the multiplication with m such that cos [marg(— Zj)\ 
can become negative. Third, with decreasing 5, the zeros 
Zj ( 5 ) are shifted towards the positive real direction in the 
complex plane. Again, deviations of arg(— Zj) from 0 be¬ 
come larger. Moreover, Eq. (2) is dominated by the zeros 
with the largest \zj(s)\~ rn . Thus, by varying 5 , different 
zeros can be brought into focus. 


III. EXAMPLES 

In the following, we illustrate our findings on a single- 
level quantum dot tunnel coupled to electronic reservoirs 


either in or out of equilibrium (cf. Fig. 1). A large 
Coulomb interaction ensures that double occupation of 
the dot is forbidden. For weak tunneling (sequential tun¬ 
neling) , electron transfer through the dot is described by 
the master equation [see Fig. 1 (c)] 


= - (7to + 740 ) Ppfit) + 7ot-fjv_i(i) + 704 .-Fat-i(*) > 


PUt) — 7fo P]sr(t) 7ot Pn(P) > 

PnP) = 74-0P}!r(t) ~ 7o±Pjv(t), 

(4) 

for the probability P$(t) that N electrons have left the 
quantum dot in time t and the dot is in a state x = 0, t? i 
(electrons entering the dot are not counted, i.e., N is not 
the net charge transfer between the dot and some lead). 
The transition rates from state x' t° X are given 
by Fermi’s golden rule. Making use of the ^-transform 
P£(t) := zN P n^), we obtain 


P*W=W*P*(t), (5) 

with the vector P z = (P^, Pj, P%) T and the matrix 


/-7to - 7to ^7ot z loi \ 

= 7to ~7ot 0 • ( 6 ) 

\ 7t0 0 “70J./ 

The solution of Eq. (5) is P z (t) = exp (W z t) P(0), 
where P(0) is the initial probability distribution. Assum¬ 
ing that electron counting starts when the system has 
reached its steady state, P(0) is given by the station¬ 
ary probability distribution, determined by WiP(0) = 
0 and e T • P(0) = 1, with e T = (1,1,1). The 
moment-generating function M s (z,t) = p z+s (t)/p s (t), 
with P z (t) = e T • P z (t), can be expressed in terms of 
the three eigenvalues A j jZ (j = 1,2,3) of the matrix 
W z = J2j ^j,z T j,z ® 1 J jZ by making use of the decom¬ 
position into left and right eigenvectors \j jZ and Tj :Z with 



FIG. 1. (Color online) (a) Equilibrium scenario: A single- 
level quantum dot subject to a Zeeman field is tunnel cou¬ 
pled to one normal lead, (b) Nonequilibrium scenario: A 
single-level quantum dot is tunnel coupled to two ferromag¬ 
netic leads with finite bias voltage, (c) Sketch of the states 
and transition rates. 
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FIG. 2. (Color online) Zeros Zj for the equilibrium scenario 
with A = k B T /2, e = —A, and times t = 5/r, 8/r, 15/r. For 
times t < t m in ~ 6.83/r, the zeros remain on the real axis. 
For t > tmin, zeros leave the real axis and interactions can be 
detected by generalized factorial cumulants. 

normalization \J Z • = 5jj> . We find 

3 

Ms(z,t) = 52 °j,z+s exp (A j, z+s t) /p a (t), (7) 

j = 1 

where we defined the amplitudes Cj, z := 

(e T • r(1 j^ z • P(0)). The generalized factorial cu¬ 

mulants are then evaluated by performing derivatives of 
In with respect to z at z = 0. 


A. Equilibrium scenario 

As a first example, we consider an equilibrium sce¬ 
nario with a quantum dot coupled to one normal lead 
[see Fig. 1 (a)]. Spin degeneracy of the dot level is lifted 
by a Zeeman field, = £± A/2, where the quantum-dot 
level 5 measured relative to the Fermi energy of the lead 
may be tuned by a gate voltage. The positive (negative) 
sign corresponds to a =t (|). Fermi’s golden rule yields 
7 cr 0 = r f(e a ) and 7 0cr = T[l-f(e a )] with the Fermi 
function f a = [l + ex.p( K £ a /kBT)]~ 1 and tunnel-coupling 
strength T. Weak tunneling corresponds to hT <C k&T. 


FIG. 3. (Color online) (Generalized) factorial cumulant C™ 
normalized by C\ — C 1 = ( N) > 0 for the equilibrium sce¬ 
nario as a function of (a) time t and (b) dot-level energy 
s. The parameters are A = k^T/2 and (a) e — —A or (b) 
t = 100/r. Negative values of (— 1 ) rn ~ 1 C™{t) indicate the 
presence of interactions. Interactions can be detected (a) for 
times larger than t m i n ~ 6.83/r and (b) for level positions 
£ < k B T\n2. 


For a vanishing Zeeman field, A = 0, the model can 
be mapped onto a noninteracting one in which a single, 
spinless quantum-dot level (x = 0,1) is filled with rate 
7io : = 7to + 7fo and emptied with rate 701 := 7ot = 704 , • 
Only two Cj^ z are nonvanishing. As a consequence, all 
zeros of M s (z,t) lie on the real axis and Eq. (3) holds. 
For a finite Zeeman field, A / 0, the zeros of M s (z,t) 
remain on the real axis for short times. 54 However, after 
a minimal time span t m i n , which depends on the system 
parameters and can be larger or smaller than 1/T, the 
first pair of zeros moves from the real axis into the com¬ 
plex plane (cf. Fig. 2). Beyond this time, the presence 
of interactions can be detected from the full counting 
statistics, as we will detail in the following. 

In Fig. 3 (a), we show the factorial cumulants (s = 1) 
as a function of time for fixed values of ejk^T and 
A/feeT. While for the first four generalized factorial cu¬ 
mulants Eq. (3) holds, there is a sign change for higher 
orders m, which indicates the presence of interactions. 
With increasing m, the time at which the sign change 
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FIG. 4. (Color online) (a) Parameter space and (b) minimal 
time span t m i n of the system in equilibrium, (a) To the left of 
the colored lines, the sign criterion for the respective C™(t) 
is violated at some time t. To the right of the dashed line, 
Ais(z,t) has only real zeros such that the sign criterion can¬ 
not be violated, (b) Minimal time span £ m in increases with 
increasing e and decreasing |A| and diverges at the dashed 
line in (a) and for A = 0. 


occurs decreases and approaches t m i n . However, when 
considering generalized factorial cumulants, we observe 
that already C B (t) violates the sign criterion Eq. (3). In 
Fig. 3 (b), the gate-voltage dependence of the factorial 
cumulants is depicted. We find that the criterion Eq. (3) 
is violated for low-lying level energies 5, indicating that 
interactions are more important in the regime when both 
spin states in the quantum dot have a finite occupation 
probability. Again, for generalized factorial cumulants 
with 8 = 0, the possibility to detect the presence of in¬ 
teractions is dramatically enhanced. 

Figure 4 (a) illustrates the possibility to detect in¬ 
teractions via different generalized factorial cumulants. 
To the right of the dashed line, given by = In 2 + 
(^^) 2 /72 + zeros °f A4 s (z,t) remain on 

the real axis [see Fig. 5 (a)] and t m i n is infinite [see 
Fig. 4 (b)]. Therefore, a violation of the sign criterion 
Eq. (3) and thus a detection of interactions is only pos¬ 
sible for e < &BTln2. The second-order factorial cu- 
mulant, 8 = 1 and m = 2, violates Eq. (3) only for 
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FIG. 5. (Color online) Zeros Zj for the equilibrium scenario 
with A = k B T /2, t = 100/r, and s/k B T — —1.5, 0.0,1.0. For 
e/k B T < In 2, zeros leave the real axis after some time and in¬ 
teractions can be detected by generalized factorial cumulants. 

rather large values of the Zeeman energy and a low-lying 
quantum-dot level. With increasing m, the region in 
which interactions can be detected is gradually increased 
to lower Zeeman splitting |A| and larger level positions 6. 
For generalized factorial cumulants with 8 = 0, already 
the second order, m = 2, covers a much larger region of 
violation of Eq. (3). This clearly demonstrates the en¬ 
hanced sensitivity of generalized factorial cumulants to 
the presence of interactions. 


B. Nonequilibrium scenario 

As a second example [see Fig. 1 (b)], we consider a 
spin-degenerate quantum dot, A = 0, coupled to two fer¬ 
romagnetic leads with parallel magnetizations and tun¬ 
nel couplings T r , r = L, R. A finite bias voltage V ap¬ 
plied symmetrically between the two ferromagnets gives 
rise to a nonequilibrium scenario. Each ferromagnet is 
characterized by its spin polarization p r ranging from 
p r = 0 for a normal metal to p r = 1 for a half- 
metallic ferromagnet with majority spins only. In the 
limit | e db eV/2| k^T, T, transport through the quan¬ 
tum dot is supported by unidirectional sequential tun- 
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FIG. 6. (Color online) (a) Parameter space and (b) minimal 
time span t m i n of the system out of equilibrium. The polar¬ 
izations are pl = p R = p, however, a different polarization 
Pl (even pl = 0) does not change the results qualitatively, 
(a) Above the colored lines, the sign criterion for the respec¬ 
tive C™(t) is violated at some time t. Below the dashed line, 
Ais(z,t) has only real zeros such that the sign criterion can¬ 
not be violated, (b) For a > —1/3, the minimal time t m i n 
increases with decreasing p and decreasing a and diverges at 
p = 0; for a < —1/3, the minimal time t m in diverges both at 
the dashed and the dotted-dashed line in (a). 


neling with rates 7 t0 = (1 + p L )r L , 7-j.o = (1 -pOU, 
7 of = (1 +Pr)Fr, and 704. = (1 - p R )r R . We denote 
the asymmetry of tunnel couplings to source and drain 

bya = (r L -r R )/(r L + r R ). 

In Fig. 6 (a), we demonstrate the possibility to de¬ 
tect the presence of interactions in this nonequilib¬ 
rium scenario. Below the dashed line, given by a = 
(3 p\ + 4 plPr + l)/(3p R - 3), all zeros of M s (z,t) re- 
main on the real axis [see Fig. 7 (c)], i.e., neither fac¬ 
torial nor generalized factorial cumulants indicate the 
presence of interactions. The same is true for the triv¬ 
ial case Pr = 0 and the dotted-dashed line, given by 
a = (Pr + 2p L p R + 1 )/(Pr~ ?PlPr ~ 3), where the sys- 
tern can be mapped onto a noninteracting Hamiltonian, 
described by a two-state model. Of course, t m j n diverges 
at these system parameters [see Fig. 6 (b)]. Above the 
dotted-dashed line, all the zeros move into the complex 
plane [see Fig. 7 (a)], including the rightmost zeros which 
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FIG. 7. (Color online) Zeros Zj for the nonequilibrium 
scenario discussed in Fig. 6 with pl = Pr = 0.35, t — 
50/(r L +r R ),anda = (r L -r R )/(r L +r R ) = -o.8,-o.5,o.o. 

For a = —0.8, all zeros remain on and for a — 0.0 all zeros 
leave the real axis. For a = —0.5, the rightmost zeros remain 
on the real axis, but zeros further to the left move into the 
complex plane. While factorial cumulants are insensitive to 
interactions in this case, generalized factorial cumulants can 
detect their presence. 


dominate the behavior of the factorial cumulants, 5 = 1. 
For large values of the spin polarization p = Pl = Pr 
or the asymmetry a of the tunnel couplings, Eq. (3) is 
already violated for the second-order factorial cumulant 
(m = 2). With decreasing p or a, higher orders m of the 
factorial cumulants are needed to detect interaction. 

An interesting regime resides between the dashed and 
dotted-dashed line. Here, the rightmost zeros remain on 
the real axis, but zeros further to the left move into the 
complex plane [see Fig. 7 (b)]. As a consequence, a viola¬ 
tion of Eq. (3) cannot be found for any ordinary (s = 1) 
factorial cumulant. This limitation is overcome by gen¬ 
eralized factorial cumulants. By an appropriate choice of 
5, they can probe the position of any zero of M s (z,t) 
in the complex plane and thus detect interactions where 
ordinary factorial cumulants fail. Furthermore, as in the 
equilibrium case, for a given order m, generalized facto¬ 
rial cumulants can detect interactions in a larger area of 
parameter space. 
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IV. CONCLUSIONS 

In summary, we proposed generalized factorial cumu- 
lants C™(t) of the full counting statistics of electron 
transport through nanostructures as a sensitive and ver¬ 
satile tool to detect the presence of interactions via the 
violation of Eq. (3). We found that generalized factorial 
cumulants are superior to ordinary ones: With decreas¬ 
ing s, interaction effects already show up in lower order 
m and at earlier times. This may be crucial for over¬ 
coming experimental limitations. Furthermore, there are 
regimes in which general factorial cumulants can detect 
interactions while ordinary ones completely fail. We il¬ 
lustrated our theoretical findings with two examples of 
a quantum dot tunnel coupled to electronic reservoirs. 
Importantly, we demonstrated that, already in a simple 


equilibrium case of experimental relevance, interactions 
can be detected via generalized factorial cumulants. Fi¬ 
nally, we emphasize the general validity of the criterion 
Eq. (3). It is not restricted to any specific type of in¬ 
teraction or transport regime and also covers multilevel 
and multichannel setups. It also applies to systems where 
coherences described by off-diagonal elements of the re¬ 
duced density matrix have to be taken into account. 
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